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Abstract 

Background: In the presence of drought and other desiccating stresses, plants synthesize and redistribute the 
phytohormone abscisic acid (ABA). ABA promotes plant water conservation by acting on specialized cells in the 
leaf epidermis, guard cells, which border and regulate the apertures of stomatal pores through which 
transpirational water loss occurs. Following ABA exposure, solute uptake into guard cells is rapidly inhibited and 
solute loss is promoted, resulting in inhibition of stomatal opening and promotion of stomatal closure, with 
consequent plant water conservation. There is a wealth of information on the guard cell signaling mechanisms 
underlying these rapid ABA responses. To investigate ABA regulation of gene expression in guard cells in a 
systematic genome-wide manner, we analyzed data from global transcriptomes of guard cells generated with 
Affymetrix ATH1 microarrays, and compared these results to ABA regulation of gene expression in leaves and other 
tissues. 

Results: The 1 173 ABA-regulated genes of guard cells identified by our study share significant overlap with ABA- 
regulated genes of other tissues, and are associated with well-defined ABA-related promoter motifs such as ABREs 
and DREs. However, we also computationally identified a unique as-acting motif, GTCGG, associated with ABA- 
induction of gene expression specifically in guard cells. In addition, approximately 300 genes showing ABA- 
regulation unique to this cell type were newly uncovered by our study. Within the ABA-regulated gene set of 
guard cells, we found that many of the genes known to encode ion transporters associated with stomatal opening 
are down-regulated by ABA, providing one mechanism for long-term maintenance of stomatal closure during 
drought. We also found examples of both negative and positive feedback in the transcriptional regulation by ABA 
of known ABA-signaling genes, particularly with regard to the PYR/PYL/RCAR class of soluble ABA receptors and 
their downstream targets, the type 2C protein phosphatases. Our data also provide evidence for cross-talk at the 
transcriptional level between ABA and another hormonal inhibitor of stomatal opening, methyl jasmonate. 

Conclusions: Our results engender new insights into the basic cell biology of guard cells, reveal common and 
unique elements of ABA-regulation of gene expression in guard cells, and set the stage for targeted 
biotechnological manipulations to improve plant water use efficiency. 




Genomics 



Background 

Drought stress imposes one of the greatest limitations to 
crop growth and yield worldwide [1]. Limitations in 
fresh water availability are predicted to become an 
increasing problem due to industrialization and continu- 
ing global climate change [2-4]. One of the central 
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mechanisms in plant drought tolerance is a reduction in 
stomatal apertures at the leaf surface, with consequent 
increase in water use efficiency [5-7]. Abscisic acid 
(ABA) is a major plant hormone which inhibits growth 
and promotes tolerance of abiotic stresses such as 
drought, salinity and cold [8-10]. Under drought condi- 
tions, plants synthesize and redistribute ABA which trig- 
gers cellular responses in guard cells, specialized cells 
that flank the stomatal pores, leading to inhibition of 
stomatal opening and promotion of stomatal closure 
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and thereby reducing plant water loss. Thus, an 
improved understanding of guard cell responses to ABA 
is relevant for development of cultivars with improved 
productivity under drought conditions. 

Guard cells also have become an advanced model sys- 
tem for understanding plant signal transduction [11-13]. 
Much of the research on guard cell responses has 
focused on the signaling cascades underlying rapid 
responses to ABA as well as to other inputs such as 
light and C0 2 [14]. However, in addition to these post- 
transcriptional responses, guard cells also respond to 
ABA at the level of the transcriptome. For example, Li 
et al. showed in Vicia faba guard cell protoplasts that 
these cells respond with an increase in dehydrin tran- 
script levels after just 15 min. of ABA exposure [15], 
suggesting that some of the cellular signaling processes 
that are observed could be influenced by transcriptional 
events. Leonhardt et al. published an expression analysis 
of the ABA-related transcriptome of guard cell proto- 
plasts [16], utilizing the Affymetrix AG GeneChip, 
which represents -8200 Arabidopsis genes. In those 
experiments, whole plants were treated with 100 uM 
ABA for four hours, and transcriptional inhibitors were 
applied during the protoplasting process, with the idea 
of stabilizing transcript levels. The authors demonstrated 
ABA-regulation of approximately 150 genes in guard 
cell protoplasts. To extend the investigation of gene 
expression in guard cells, Yang et al. isolated guard cell 
promoter candidates based on 4 Arabidopsis guard cell 
microarrays from 4 different treatments using whole- 
genome Affymetrix ATH1 chips [17], which profile 
-24,000 genes. 

To date, however, there has been no systematic repli- 
cated study that investigates ABA regulation of gene 
expression in guard cells in a genome-wide manner. A 
comprehensive analysis of the global ABA-related tran- 
scriptome of guard cells and its relationship to known 
guard-cell signaling cascades promises to engender new 
insights into the basic cell biology of guard cells. In 
addition, the wealth of microarray data regarding ABA- 
regulated transcriptomes of other tissues, including 
whole plants [18,19], seedlings [20-26], and seeds [27], 
provides an opportunity to assess the extent to which 
the specialized guard cells employ unique vs. common 
promoter motifs and unique vs. common gene targets in 
ABA regulation of gene expression. Given growing evi- 
dence of hormonal cross-talk in physiological responses 
[10], another interesting question is whether or not 
ABA-regulated transcripts of guard cells are regulated 
by other hormones in other plant tissues. Finally, 
manipulation of transcriptional regulation of key genes 
is one of the most common approaches utilized in 
applied research aimed at improving plant stress toler- 
ance [28-30], so genome-wide analysis of the 



Arabidopsis guard cell transcriptome may yield candi- 
date genes for translational research. In this study, we 
performed a systematic microarray analysis of gene 
expression profile data from global transcriptomes of 
guard cells generated with Affymetrix ATH1 chips to 
explore the above topics. 

Results 

Identification of ABA-responsive genes in guard cells and 
leaves 

We and others [16] have observed that when guard cell 
protoplasts were used as the source of RNA, stress- 
induced genes could be induced by protoplasting, even 
without ABA treatment. Therefore, we employed a pro- 
tocol [31] whereby isolated epidermal peels with guard 
cells as the only intact cell type (Figure 1A and Figure 
IB) were used as the source of guard cell RNA for 
microarray analysis. 

For each tissue by treatment combination (guard cell 
vs. leaf; ABA vs. ethanol as solvent control), three inde- 
pendent biological samples were assessed using Affyme- 
trix whole genome ATH1 chips [31]. For the present 
study, to reliably identify ABA-responsive genes, we 
adopted two different methods - a Boolean method [31] 
and a standard linear model method [32] - and focused 
on the set of ABA-regulated genes jointly confirmed by 
both approaches (see Methods for details). Through this 
overlapping analysis we identified 1173 ABA-responsive 
genes in Col (Arabidopsis Col-0) guard cells (696 ABA- 
induced genes and 477 ABA-repressed genes) and 1037 
ABA-responsive genes in Col leaves (596 ABA-induced 
genes and 441 ABA-repressed genes). The ABA-regu- 
lated genes in guard cells and leaves are listed in Addi- 
tional file 1 and Additional file 2, respectively. As shown 
in Figure 1C and ID, among all the ABA-induced genes, 
only 177 genes are ABA-induced in both guard cells and 
leaves, while the rest show differential regulation with 
respect to tissue type. The overlap of guard cell ABA- 
repressed genes and leaf ABA-repressed genes is even 
smaller with only 90 genes common to both guard cells 
and leaves. There are only 4 genes that are oppositely 
regulated by ABA in guard cells vs. leaves: AT5G20840, 
EXL3 (AT5G51550), and CIL (AT4G25990) are induced 
by ABA in guard cells but repressed by ABA in leaves, 
and KCSS (AT1G25450) is induced by ABA in leaves but 
repressed by ABA in guard cells. 

Validation of ABA-responsive genes from microarray 
analysis by Q-PCR 

To evaluate the reliability of ABA-regulation derived 
from the microarray data and validate identified ABA- 
responsive genes, we performed Q-PCR experiments on 
3 genes identified as not responsive to ABA and 14 genes 
that showed differential expression in response to ABA 
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Figure 1 Images of guard cells and overview of guard cell ABA-regulated genes and leaf ABA-regulated genes. (A) Bright-field image: 
Epidermal peel with intact guard cells after partial cell wall digestion. (B) Fluorescent image: Vital staining of epidermal peels with FDA reveals 
that guard cells are the only intact cell type. (C) Overlap of guard cell ABA-induced genes and leaf ABA-induced genes. (D) Overlap of guard cell 
ABA-repressed genes and leaf ABA-repressed genes. 



in either guard cells, leaves, or both tissues by our micro- 
array analysis (Figure 2). The three negative control 
genes were Actin {Actin 2: AT3G18780 and Actin 8: 
AT1G49240), a /3-tubulin like protein TUB9 
(AT4G20890), and a putative nematode resistant protein 
Arabidopsis ortholog of sugar beet HS1 PRO-1 2 
{HSPR02, AT2G40000). Among the ABA-responsive 
genes, we chose two genes that exhibited ABA-upregula- 
tion in both guard cells and leaves in our microarray ana- 
lysis: ABI2 (AT5G57050), which encodes one of several 
PP2C-type protein phosphatases that negatively regulate 
many aspects of ABA signaling, including ABA-induced 
stomatal closure [33-35], and the non-specific lipid trans- 
fer protein 3 (LTP3, AT5G59320), which is a drought- 
responsive gene [36] . In addition, we analyzed two genes 
that showed ABA-repression in both guard cell and leaf 
microarrays: a heat shock protein (AT4G21870) and a 
somatic embryogenesis receptor-like kinase 2 (SERK2, 
AT1G34210). For the category of ABA-induction in 
guard cell but not leaf microarrays, five genes were cho- 
sen for Q-PCR evaluation: LTI78/RD29A (AT5G52310), 
identified as a cold and salt stress-induced gene in a leaf 



transcription study [37,38], two genes encoding phospho- 
lipases known to be involved in ABA signaling, namely a 
phospholipase C gene, PI-PLC (AT3G55940) and a PLD8 
homolog (AT4G35790), as well as GBF2 (AT4G01120), 
and a protein kinase (AT1G77280). For genes that 
showed ABA-repression in guard cell but not leaf micro- 
arrays, we performed Q-PCR assessment of a vacuolar 
anion channel gene ATCLC-A (AT5G40890) and the ion 
channel gene KAT1 (AT5G46240), which mediates K + 
uptake during stomatal opening. For genes that were 
ABA-responsive only in the leaf microarrays, we tested 
ECERIFERUM 1 {CER1, AT1G02205) and Ca 2+ -binding 
protein 1 (ATCP1, AT5G49480) for ABA-induced genes 
and WALL-ASSOCIATED KINASE 2 ( WAK2, 
AT1G21270) for ABA-repressed genes. 

Figure 2 summarizes the Q-PCR analysis. For 15 out 
of 17 genes, the Q-PCR results match the microarray 
data. The two exceptions are LTI78/RD29A and PI-PLC, 
for which ABA-upregulation is observed in guard cells 
but not in leaves from the microarray analysis, while the 
Q-PCR results show ABA-upregulation in both guard 
cells and leaves. This 94.1% (= 32/34 comparisons) 
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Figure 2 Validation of a set of ABA responsive genes identified from microarray analysis by Q-PCR GC: guard cells. LF: leaves (A) 
Expression fold changes of the reference gene Actin and two control genes in response to ABA. (B) Expression fold changes of genes that are 
up-regulated by ABA in both leaf and guard cell microarrays. (C) Expression fold changes of genes that are down-regulated by ABA in both leaf 
and guard cell microarrays. (D) Expression fold changes of genes that are up-regulated by ABA in guard cell microarrays only. (E) Expression fold 
changes of genes that are down-regulated by ABA in guard cell microarrays only. (F) Expression fold changes of genes that are up-regulated by 
ABA {ATCP1 and CER1) and down-regulated by ABA (WAK2) in leaf microarrays only. See text for specific gene locus identifiers. 



agreement between the microarray analysis and the Q- 
PCR experiments is comparable or better than that 
observed in other microarray studies [39]. The discre- 
pancy most probably stems from the different techni- 
ques employed to measure gene expression. Several 



studies have noted that microarray analysis is less sensi- 
tive than Q-PCR and can underestimate the true extent 
of transcript accumulation [40,41], so observing this 
result for a minority of the genes (2 out of 17) is not 
surprising. 
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c/s-acting regulatory elements in ABA-responsive genes 

Several c/s-regulatory elements responsible for mediating 
ABA-, drought-, and cold-induced gene expression have 
been identified by biochemical and molecular techniques 
[42]. We were interested to determine whether these 
known motifs would be enriched in our gene sets, 
which would support our identification of these genes as 
ABA-regulated. In addition, we were interested to assess 
whether ABA-regulated genes of guard cells and leaves 
would exhibit differential usage of these promoter 
motifs, or differential positioning of the motifs within 
the promoter region. 

The G-box-containing ABA-responsive elements 
(ABREs) and coupling element 3 (CE3) are classical ex- 
acting elements in ABA-responsive gene expression 
[43-46]. In addition, the core motif of the dehydration- 
responsive element/C-repeat (DRE/CRT) is a c/s-acting 
element found in many drought-, high-salt- and cold- 
responsive genes in Arabidopsis and rice, and provides 
the binding site for DRE-related transcription factors 
such as DRE-Binding proteins (DREBs) [46,47]. MYB 
and MYC binding motifs are also found in the promo- 
ters of drought-inducible genes. Proteins of the MYB 
and MYC families function in mediating drought- and 
ABA-regulated gene expression [33,48] and some mem- 
bers of these families such as AtMYC2 and AtMYB2 are 
synthesized following accumulation of endogenous ABA 
[49]. Other related c/s-acting elements include the bind- 
ing site of C-repeat factors (CBFs) and the low tempera- 
ture responsive element (LTRE) [46]. 

We performed a statistical analysis (see Methods) of the 
promoter sequences (defined as the 1000-bp upstream 
regions) of ABA-responsive genes in guard cells and 
leaves, to check whether the known motifs would be 
enriched in our gene sets. The results of this analysis are 
shown in Table 1. We can see that ABREs are significantly 
enriched in guard cell ABA-induced genes. In particular, 
more than 50% of ABA-induced genes in guard cells have 
ABREs in their promoter regions, supporting the reliability 



of our microarray methods and predictions. DREs and 
LTREs are also highly enriched in guard cell ABA-induced 
genes. ABREs and the cold-related c/s-acting element 
LTRE, but not DRE/CRT elements are also highly 
enriched in leaf ABA-induced genes. A motif over-repre- 
sented in light-induced promoters, GCCAC (SORLIP1AT 
in the PLACE database [46]), also is found to be highly 
enriched in guard cell ABA-regulated genes but not in 
leaves. Compared with other motifs, CE3 has only a few 
matches in the ABA-regulated genes of either sample type. 
This is consistent with the findings that CE3 is well-repre- 
sented in rice but almost absent in Arabidopsis [50]. None 
of the as-acting elements enriched in ABA-induced genes 
is significantly enriched in ABA-repressed genes of guard 
cells or leaves, except the MYC binding site and LTRE. 
The mechanisms of ABA-repression are less well charac- 
terized than those for induction [33] and motifs specific 
for ABA-repressed genes may exist but have not been 
revealed yet (see next section). 

We also examined the positional preference of the c/s- 
acting regulatory elements ABRE, DRE, CBF and LTRE 
in a 1000-bp region upstream of identified ABA-regu- 
lated genes of guard cells and leaves (Figure 3). There is 
no significant difference between ABRE locations in 
guard cell vs. leaf genes (Figure 3A); the same is true for 
the c/s-acting elements DRE, CBF and LTRE locations in 
guard cell genes and leaf genes (not shown). The posi- 
tional distribution of these four elements upstream of 
guard cell ABA-induced genes is shown in Figure 3B. 
ABREs are more likely than the other three motifs to be 
located in regions close to transcription start sites (TSS), 
with around 28% and 20% ABREs located in the -100— 
50-bp region upstream of the TSS in guard cell and leaf 
ABA-induced genes, respectively. 

Identification of new promoter motifs 

In guard cells, 39 of the ABA-induced genes have none 
of the known motifs of Table 1 in their promoter 
regions, and this is also true for 79 of the genes induced 



Table 1 Enrichment of known c/s-acting regulatory elements in 1000-bp region upstream of ABA-regulated genes 



Motif name 


Motif sequence 


All genes 


GC„up (696) 


GC„down (477) 


LF_up (596) 


LF_down (441) 


Hits 


Hits 


P-value 


Hits 


P-value 


Hits 


P-value 


Hits 


P-value 


ABRE 


(C/G/T)ACGTG(G/T)(A/C) 


3178 


382 


7.1E-149 


49 


0.99 


188 


2.3E-29 


57 


0.71 


CE3 


ACGCGTGTC 


24 


7 


3.5E-7 


0 


040 


6 


1.9E-6 


I 


0.08 


DRE/CRT 


(A/QCCGAC 


2748 


137 


2.0E-9 


56 


0.56 


85 


0.04 


58 


0.22 


CBF binding 


(A/G)(C/T)CGAC 


6236 


284 


3.6E-15 


131 


0.47 


215 


1.2E-6 


133 


0.09 


LTRE 


CCGAC 


5062 


229 


2.0E-1 1 


106 


0.48 


162 


1.8E-3 


119 


7.7E-3 


MYB binding 


{C/T)AAC(G/T)G 


9108 


313 


3.2E-3 


213 


0.02 


246 


0.25 


195 


0.03 


MYC binding 


CACATG 


4909 


180 


2.7E-3 


138 


5.5E-5 


165 


1.5E-4 


118 


3.8E-3 


SORLIP1 AT 


GCCAC 


6609 


321 


1.4E-22 


152 


0.08 


189 


0.07 


136 


0.19 



'All genes' means all genes represented in the Affymetrix ATH1 chips. "GC_up" and "GC_down" represent ABA-induced and ABA-repressed genes in guard cells, 
respectively. "LF_up" and "LF_down" mean ABA-induced and ABA-repressed genes in leaves, respectively. The numbers in brackets in the header denote the 
total number of ABA-regulated genes in the corresponding category. P-values smaller than 0.01 are marked in bold. 
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by ABA in leaves. In addition, not all genes containing a 
known ABA-related promoter motif show ABA-regula- 
tion in both tissue types. Accordingly, we sought to 
identify possible new motifs that may be involved in 
ABA-regulated gene expression, particularly of a cell- 
specific nature. We were also interested in the possibi- 
lity of identifying repression-specific motifs. We identi- 
fied all 5-10mer motifs highly enriched in the 1000-bp 
sequences upstream of ABA-regulated genes in guard 
cells and leaves (see Methods), and also compared their 
enrichment in previously published ABA-regulated gene 
sets. The statistically significant 5-mer motifs are sum- 
marized in Table 2. These motifs include not only the 
core of G-box-containing ABREs (ACGTG), subse- 
quences of ABREs, the LTRE (CCGAC), and the light- 
induced motif (GCCAC), as expected from the above 
analysis, but also some new, previously unidentified, 
motifs. The results for 5-10mer motifs also show that 
there are many more significant motifs in ABA-induced 
genes than in ABA-repressed genes. In addition, a 
greater number of distinct motifs were found in guard 
cell ABA-induced genes than in leaf ABA-induced 
genes. 

To check whether newly identified 5-mer motifs sig- 
nificantly enriched in our gene sets are guard cell- or 
leaf-specific, we merged previously published ABA-regu- 
lated gene sets in Arabidopsis from [18-27] and com- 
pared the enrichment of these motifs. Note that none of 
these studies uses guard cells or rosette leaves as plant 
materials. We calculated the enrichment P-values of 
these motifs in the merged ABA-regulated gene set 
using all genes in the whole genome as the background, 
and also calculated the enrichment of these motifs in 
our gene sets using the merged ABA-regulated gene set 
as the background. The results are shown in Table SI 



and Table S2 in Additional file 3, respectively. Guard 
cell-specific or leaf-specific motifs are designated as 
those that are not significantly enriched in the above 
merged ABA-regulated gene set, but are still significant 
in our guard cell or leaf ABA-regulated gene set when 
using the merged ABA-regulated gene set as the back- 
ground. We find that GTCGG, which is not a subse- 
quence of any known motifs, is significantly enriched in 
guard cell ABA-induced genes using either all genes in 
the whole genome or the merged ABA-regulated gene 
set as the background. By contrast, the enrichment of 
GTCGG in the merged ABA-regulated gene set is not 
significant (larger than 1.0E-03). We accordingly identify 
GTCGG as a potential guard cell-specific motif involved 
in ABA-induced gene expression. Similarly, TGCAA is a 
potential leaf-specific motif involved in ABA-repressed 
gene expression. 

Although only a few motifs emerge as significant from 
the promoters of ABA-repressed genes, in our analysis 
we find a novel strong candidate for a repression-speci- 
fic motif, CAAGTTG, which is enriched in both guard 
cell ABA-repressed genes and leaf ABA-repressed genes 
(P = 6.0E-09 and 1.6E-08, respectively). This motif is 
very similar to the motif CAACTTG identified in [51], 
but in our guard cell and leaf gene sets, the enrichment 
of CAAGTTG is stronger than that of CAACTTG (P = 
6.6E-05 and 1.7E-04, respectively). Known ABA- 
repressed genes KAT2 and MYB60 [14,16] contain the 
CAAGTTG motif. To examine whether this motif is 
also enriched in ABA-repressed genes in other tissues/ 
cell types, we collected eight ABA-repressed gene sets 
from [19-26], and compared this motifs enrichment in 
these ABA-repressed gene sets (studies that report a 
very small number of ABA-repressed genes are not 
included here). As shown in Figure 4, CAAGTTG is 
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Table 2 5-mer motifs significantly enriched in 1000-bp regions upstream of guard cell or leaf ABA-regulated genes 

5-mer motifs ABA-induced genes 5-mer motifs ABA-repressed genes 



GC (P< 10" 10 ) LF(P < 10" 10 ) GC (P< 10" 4 ) 3 LF(P < 10" 4 ) a 



ACGTG 


7.5E-90 


4.2E-32 


CCACT 


4.3E-07 


2.6E-04 


CACGT 


1 .2E-84 


1 .3E-33 


CCAAC 


2.1 E-06 


9.6E-04 


CCACG 


5.7E-56 


4.1 E-1 1 


CAACT 


3.0E-05 


1 4E-04 


ACACG 


1.5E-51 


5.8E-21 


CACAT 


7.2E-05 


1 .1 E-03 


CGTGT 


1 .7E-47 


1 .1 E-22 


GGTCC 


1 .3E-04 


7.7E-07 


CGTGG 


5.6E-36 


2.0E-08 


TGCAA 


2.5E-03 


1 .1 E-05 


GTGTC 


1 .7E-24 


2.5E-03 


GTCCC 


7.5E-04 


2.8E-05 


GCCAC 


8.6E-23 


>0.01 


GACCA 


>0.01 


5.3 E-05 


GTGGC 


1.3E-16 


5.9E-05 








GACAC 


5.9E-15 


>0.01 








GACGT 


7.3E-13 


3.1E-07 








CACGC 


1.2E-12 


1.2E-05 








TACGT 


5.5E-12 


6.6E-12 








GTCGG 


9.9E-12 


6.8E-03 








CCGAC 


1.5E-11 


1.5E-03 








ACGTA 


1.7E-11 


1.3E-13 








TCCAC 


2.9E-1 1 


2.2E-05 








TGTCG 


8.9E-1 1 


5.8E-03 









Motifs marked in bold are subsequences of ABRE considered in Table 1. Underlined motifs exhibit tissue specificity (GC vs LF) based on our analysis (Additional 
file 3). For convenient comparison of GC and LF, P-values falling beyond the thresholds (see Methods) are also listed. a Motifs appearing in more than 15000 loci 
are not listed. 



Arabidopsis tissues and cell types, including whole plants 
[18,19], seedlings [20-26], seeds [27], and one previous 
study reporting ABA responsive genes in guard cells and 
mesophyll cells using a partial-genome (~8K) Affymetrix 
chip [16]. 5789 ABA-induced genes and 5635 ABA- 
repressed genes are obtained by merging the ABA- 
responsive genes in these studies. We examined how 
many of the ABA-regulated genes that we identified from 
our guard cell transcriptomes can be found in this 
merged ABA-regulated gene set. The overlap of our 
ABA-regulated genes with the previously published 
ABA-regulated genes is shown in Figure 5. We find that 
69.4% of the ABA-induced genes and 54.5% ABA- 
repressed genes in guard cells identified in this work are 
found in previous ABA transcriptome studies. In leaves, 
86.6% of ABA-induced genes and 58.0% ABA-repressed 
genes are found in previous ABA transcriptome studies. 
Such significant overlap supports the reliability of our 
microarray analysis. 

In addition to the overall comparison, we also con- 
ducted a pairwise comparison of the 11 ABA transcrip- 
tome studies mentioned above and our study, which 
describe a total of 14 different ABA-regulated transcrip- 
tomes (experiments) (Leonhardt et al. [16] contribute 
two experiments, and we contribute two as well). These 
transcriptome experiments have been done under var- 
ious experimental conditions, including different ABA 
doses and treatment durations, tissue/cell types, micro- 
array platforms, plant developmental stages, and 



highly enriched in the upstream sequences of most of 
the eight ABA-repressed gene sets, and the combined 
motif CAA(G/C)TTG is consistently enriched in all 
ABA-repressed gene sets, which further supports its 
involvement in ABA-mediated gene repression. 

Comparison with other ABA transcriptome studies 

Many transcriptome studies have been conducted on 
ABA regulation of gene expression in different 




E1 E2 E3 E4 E5 E6 E7 E8 E9 E10 
ABA-repressed gene sets 

Figure 4 Enrichment of the CAAGTTG motif, the CAACTTG 
motif and the CAA(G/C)TTG motif in 1000-bp upstream 
sequences of the ABA-repressed gene sets from ten 
transcriptomes. El: Hoth ef at, 2002 [20]; E2: Li et al, 2006 [21]; E3: 
Xin ef al., 2005 [25]; E4: Matsui ef al, 2008 [19]; E5: Huang et al, 
2007 [26]; E6: Zeller et al., 2009 [23]; E7: Sanchez et al, 2004 [24]; E8: 
Nemhauser et al, 2006 [22]; E9: this study (guard cells); E10: this 
study (leaves). The dotted line indicates the threshold for 
significance (P < 0.01) 
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A B 

Figure 5 Overlap of ABA-regulated genes in guard cells and 
leaves with previously published ABA-regulated gene sets (A) 

ABA-induced genes. (B) ABA-repressed genes. Studies contributing 
to the merged gene set were: Hoth ef al, 2002 [20], Li ef al., 2006 
[21], Xin ef ah, 2005 [25], Matsui ef al., 2008 [19], Huang ef al., 2007 
[26], Zeller ef al, 2009 [23], Sanchez ef al, 2004 [24], Nemhauser ef 
al., 2006 [22], Leonhardt ef al, 2004 [16], Seki ef al., 2002 [18], 
Okamoto ef al., 2010 [27]. 



methods for identifying ABA-regulated genes. A sum- 
mary of experimental conditions for these ABA tran- 
scriptome studies is given in Table S3 (Additional file 
3). The pairwise overlaps are measured by percentage 
(Additional file 3) and representation factor (see Meth- 
ods), and their significance is examined by the hyper- 
geometric distribution. The results are shown in Figure 
SI (Additional file 3). Given the diversity of experimen- 
tal parameters, it is not surprising that these studies do 
not have very high pairwise overlap percentages for 
ABA-regulated genes. However, when compared with 
random overlaps, all of them are significant for ABA- 
induced genes. For ABA-repressed gene sets, all have 
significant pairwise overlaps as compared with the 
expected overlap of two random sets except for the 
ABA-repressed gene sets in [16,18,27] which have non- 
significant overlaps with several other studies. Figure 6 



shows the cumulative number of ABA-regulated genes 
that are common to different numbers of experiments. 
As expected from the pairwise overlap analysis, there 
are not many ABA-induced genes found in common as 
the number of experiments increases, and there are no 
ABA-repressed genes that are common to >11 experi- 
ments (transcriptomes). However, despite the diversity 
in tissue types, ABA treatments, and transcriptome plat- 
forms utilized in these experiments, we were able to 
find a "core set" of -50 ABA-induced genes reported in 
any 9 out of 14 experiments (Figure 6A), and a similar 
number of "core set" ABA-repressed genes reported in 
any 6 out of 14 experiments (Figure 6B). Table 3 lists 
the core set of ABA-induced genes, and Table 4 lists the 
core set of ABA-repressed genes. 

Identification of new guard cell ABA-responsive genes 

In addition to identifying core sets of ABA-induced and 
ABA-repressed genes, we also discovered a number of 
novel ABA-responsive genes in both guard cells and 
leaves. In particular, our results provide a valuable com- 
prehensive source of ABA-regulated genes in guard cells 
since the smaller Affymetrix ~8K array was used in [16] 
and therefore could not provide information on -2/3 of 
the transcriptome. For example, among 1173 genes sig- 
nificantly ABA-regulated in our guard cell transcrip- 
tome, 1122 genes (658 ABA-induced and 464 ABA- 
repressed) were not previously reported to be ABA- 
regulated by Leonhardt et al. [16]. Of these 1173 genes, 
148 ABA-induced genes and 149 ABA-repressed genes 
were also not reported as ABA-regulated in any of the 
11 previous transcriptome studies we analyzed (see Fig- 
ure 5 legend for a list of these studies). The proposed 
guard cell-specific ABA-induction motif, GTCGG, is sig- 
nificantly enriched {P < 1.3E-04) in the 148 newly 
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Table 3 67 ABA-induced genes common to > 9 experiments 


Genes 


Descriptions 


No. 


AT1G01470* A 


LEA 14 (LATE EMBRYOGENESIS ABUNDANT 14) 


14 


AT2G33380* A 


RD20 (RESPONSIVE TO DESSICATION 20); calcium ion binding 


13 


AT3G1 1410* A 


AHG3/ATPP2CA (ARABIDOPSIS THALIANA PROTEIN PHOSPHATASE 2CA) 


13 


AT1 G20440* A 


COR47 (cold regulated 47) 


12 


AT1 G52690* A 


late embryogenesis abundant protein, putative 


12 


AT1G72770* A 


HAB1 (HOMOLOGY TO ABM); protein serine/threonine phosphatase 


12 


AT2G41 1 90* A 


amino acid transporter family protein 


1 2 


AT2G47770* 


benzodiazepine receptor-related 


1 2 


AT5G06760* A 


Ig^g embryogenesis abundant group 1 domain-containing protein 


12 


AT1 G08920* 


sugar transporter, putative 




AT1 G77450* 


ANAG032 (Arabidopsis NAC domain containing protein 32) 


! ] 


AT2G15970* A 


COR413-PM1 (cold regulated 413 plasma membrane 1) 


! | 


AT2G46270* A 


GBF3 (G-BOX BINDING FACTOR 3); transcription factor 




AT2G46680 A 


ATHB-7 (ARABIDOPSIS THALIANA HOMEOBOX 7); transcription factor 




AT4G27410* 


RD26 (RESPONSIVE TO DESSICATION 26); transcription factor 


! ! 


AT5G52310* A 


COR78 (COLD REGULATED 78) 


! ! 


AT5G59370* 

r \ 1 JvJ JZ.U 


I TP3 (I IPID TRANSFFR PROTFIN 3V linid bindinn 


! } 


AT5G66400* A 


RAB18 (RESPONSIVE TO ABA 18) 


! } 


AT1 G04220* A 


KCS2 (3-KETOACYL-COA SYNTHASE 2); fatty acid elongase 


10 


AT1G05100* A 


MAPKKK18 (Mitogen-activated protein kinase kinase kinase 18) 


10 


ATI G 16850* 


unknown protein 


10 


AT1 G20450 


ERD10/LTI45 (EARLY RESPONSIVE TO DEHYDRATION 10) 


10 


ATI G49450* 


transducin fami y protein/WD-40 repeat family protein 


10 


AT1 G51 140* 


basic helix-loop-helix (bHLH) family protein 


10 


AT1 G58270* 


ZW9 


10 


AT3G1 1420 


fringe-related protein 


10 


AT3G29575* 


AFP3 (ABI FIVE BINDING PROTEIN 3) 


10 


AT3G61890 A 


ATHB-12 (ARABIDOPSIS THALIANA HOMEOBOX PROTEIN 12) 


10 


AT4G05100* 


AtMYR74 rmvh Hnm^in nrntpin 74V DMA hinHinn/fr^nsrnntinn fartnr 


10 


AT4G24130* A 


unknown protein 


10 


AT4G26080* 


ABI1 (ABA INSENSITIVE 1); calcium ion binding 


10 


AT4G30960* 


CIPK6 (CBL-INTERACTING PROTEIN KINASE 6); kinase 


10 


AT5G01520* 


zinc finger (C3HC4-type R NG finger) fami y protein 


10 


AT5G 15960 


Stress-induced protein KIN1 


10 


AT5G59220* 


protein phosphatase 2C, putative/PP2C, putative 


10 


AT1 G07430* 


protein phosphatase 2C, putative 


g 


AT1 G07720 


KCS3 (3-KETOACYL-COA SYNTHASE 3) 


9 


AT1G21790* 


similar to unnamed protein product [Vitis vinifera] (GB;CA061 872.1) 


9 


AT1G60190* 


armadillo/beta-catenin repeat family protein/U-box domain-containing protein 


9 


AT1 G62570* 


FMO GS-OX4 (FLAVIN-MONOOXYGENASE GLUCOSINOLATE S-OXYGENASE 4) 


9 


AT1G62710* 


BETA-VPE (vacuolar processing enzyme beta); cysteine-type endopeptidase 


9 


AT1 G 73480* 


hydrolase, alpha/beta fold family protein 


9 


AT1G77120* 


ADH1 (ALCOHOL DEHYDROGENASE 1) 


g 


AT2G 17840 A 


ERD7 (EARLY-RESPONSIVE TO DEHYDRATION 7) 


9 


AT2G30360* 


CIPK11 (SOS3-INTERACTING PROTEIN 4); kinase 


9 


AT2G30550* 


lipase class 3 family protein 


9 


AT2G37870 


protease inhibitor/seed storage/lipid transfer protein (LTP) family protein 


9 


AT2G39050* 


hydroxyproline-rich glycoprotein family protein 


9 


AT2G47780* 


rubber elongation factor (REF) protein-related 


9 


AT3G02480* 


ABA-responsive protein-related 


9 


AT3G17520* A 


late embryogenesis abundant domain-containing protein 


9 
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Table 3 67 ABA-induced genes common to ?>? 9 experiments (Continued) 



A 1 3bzzoUU 


protease inhibitor/seed storage/ ipid transfer protein (LTP) family protein 


9 


A I 3lj4o3 I U 


unknown protein 


y 


A I 3bDUy/U A 


LMJU (LUVV 1 tlVlr tnA 1 Unt-IINUULtU jUJ 


Q 

y 


A I jujjjUU 


ATCVDM^ /ADADinnDCIC Tl— 1 A 1 IAMA EVDAMCIM A 1 C\ 

AltArAlo (AKADlUUrblb I HALIAINA tArAIMblN Am} 


9 


A I Jbz>/U I U 


strictosidine synthase family protein 


9 


A 1 4b I /ioU 


transporter-related 


y 


A I 4uz I 44U A 


ATM/l /ATI\/1VD1 m / A D A PI HODCIC MVR I ll/C miV tmnrrrintinn fn«-t/-\r 

A I IV14/A uviyd I uz (AnABi uuk bi j iviTD-Liixt I uz), transcription ractor 


y 


AT4G23050* 


protein kinase, putative 


9 


AT4G30470* 


cinnamoyl-CoA reductase-related 


9 


AT4G33550 


lipid binding 


9 


AT4G34000*a 


ABF3/DPBF5 (ABSCISIC AGID RESPONSIVE ELEMENTS-BINDING FACTOR 3) 


9 


AT5G11110* 


ATSPS2F/SPS1 (SUGROSE PHOSPHATE SYNTHASE 1) 


9 


AT5G52300*a 


LTI65/RD29B (RESPONSIVE TO DESSICATION 29B) 


9 


AT5G57050* 


ABI2 (ABA INSENSITIVE 2); protein serine/threonine phosphatase 


9 


AT5G59310* 


LTP4 (LIPID TRANSFER PROTEIN 4); lipid binding 


9 


AT5G61820* 


similar to MtNl 9-like protein [Pisum sativum] (GB:AAU 14999.2) 


9 



The first column lists the AGI locus identifiers of the genes. *Genes that appear in our guard cell ABA-induced gene set. A Genes that appear in Leonhardt et al.'s 
guard cell ABA-induced gene set [16]. The genes that appear in our ABA-induced gene set are marked with underlines if they are not represented on the 
Affymetrix ~8K chips used in [16]. The second column gives brief TAIR descriptions for the genes. The third column indicates the number of experiments that the 
genes are found common to. 



identified ABA-induced genes of guard cells, consistent 
with the hypothesis that this motif participates in ABA- 
regulation of gene expression in guard cells. We also 
identified a smaller set of ABA-regulated transcripts in 
leaves that have not been reported in the previous ABA 
microarray studies: 62 ABA-induced genes and 113 
ABA-repressed genes. In Table 5 and Table 6 we list by 
.P-value the top 50 new ABA-induced genes and ABA- 
repressed genes in guard cells not previously reported 
by any of the 11 previous ABA microarray studies that 
we evaluated. Additional file 4 and Additional file 5 list 
all the new ABA-responsive genes in guard cells and 
leaves, respectively. 

Functional enrichment in ABA-responsive genes 

ABA is involved in a myriad of biological processes, 
including plant growth and development, and plant 
responses to environmental stresses such as cold, 
drought, and pathogens. To examine whether ABA- 
responsive genes identified in our study tend to be 
involved in such aspects of plant physiology, we used 
the GO functional analysis tool BiNGO 2.42 embedded 
in the Cytoscape project [52,53] to functionally categor- 
ize our ABA-responsive genes in guard cells and leaves. 
TAIR9 GO data (GOSlim_Plants) [54] for the whole 
Arabidopsis gene annotation were used as the reference 
set. The significance of overrepresentation and underre- 
presentation of a functional category was calculated by 
the hypergeometric cumulative distribution. Bonferroni 
familywise error rate (FWER) correction was applied for 
multiple testing. The functional categories with P-values 



less than 1.0E-03 that also have FWER -corrected P- 
values smaller than 0.05 are given in Figure 7. 

For the ABA-induced genes, the three most significant 
functional categories in both guard cells and leaves are 
response to abiotic stimulus, response to endogenous sti- 
mulus, and response to stress (Figure 7A), all consistent 
with the fact that ABA plays important roles in the reg- 
ulation of gene expression during stress responses. Lipid 
metabolic process, catalytic activity and secondary meta- 
bolic process are also enriched in both guard cell and 
leaf ABA-induced genes but the enrichment of vacuole 
and transcription factor activity is unique to guard cells. 
In addition, several functional categories are underrepre- 
sented within ABA-induced genes, including translation, 
RNA binding, protein metabolism, and nucleobase meta- 
bolism (see Discussion). 

ABA-repressed genes show enrichment of different 
functional categories than ABA-induced genes (Figure 
7B). The functional categories enriched in both guard 
cell and leaf ABA-repressed genes are plastid, transferase 
activity, plasma membrane, and nucleotide binding. 
Response to abiotic stimulus, response to stress and several 
other functional categories are highly enriched in guard 
cell ABA-repressed genes but not in leaf ABA-repressed 
genes. Two functional categories, nucleobase metabolic 
process and protein modification process, are enriched in 
leaf ABA-repressed genes but not in guard cell ABA- 
repressed genes. The functional categories protein meta- 
bolism, nucleobase metabolism and nucleotide binding 
which are underrepresented in leaf ABA-induced genes 
are overrepresented in leaf ABA-repressed genes. 
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Table 4 47 ABA-repressed genes common to > 6 experiments 


Genes 


Descriptions 


No. 


AT2G38310* 


unknown protein 


10 


AT4G21870* 


26.5 kDa class P-related heat shock protein (HSP26.5-P) 


10 


AT1 G66940* 


protein kinase-re ated 


9 


AT1 G69530* 


ATEXPA1 (ARABIDOPSIS THALIANA EXPANSIN A1) 


9 


ATI G03870* 


FLA9 


8 


AT1 G08930 


ERD6 (EARLY RESPONSE TO DEHYDRATION 6) 


8 


AT4G36670* 


mannitol transporter, putative 


8 


AT5G02760* 


protein phosphatase 2C family protein 


8 


AT5G05440* 


unknown protein 


8 


AT1G14210* 


ribonuclease T2 family protein 


7 


AT3G01860* 


unknown protein 


7 


AT3G 14840* 


eucine-rich repeat family protein/protein kinase family protein 


7 


AT3G50740 


UGT72E1 (UDP-glucosyl transferase 72E1) 


7 


AT4G 17460* 


HAT1 (homeobox-leucine zipper protein 1); DNA binding/transcription factor 


7 


AT1 G07090 


LSH6 (LIGHT SENSITIVE HYPOCOTYLS 6) 


6 


AT1G08810* A 


MYB60 (myb domain protein 60); DNA binding/transcription factor 


6 


AT1 G29430 


auxin-responsive fami y protein 




AT1G51850 


eucine-rich repeat protein kinase, putative 


6 


AT1G52190 


proton-dependent oligopeptide transport (POT) fami y protein 




AT1 G56430 


NAS4 (NICOTIANAMINE SYNTHASE 4); nicotianamine synthase 




AT1 G68840 


RAV2 (REGULATOR OF THE ATPASE OF THE VACUOLAR MEMBRANE) 


6 


AT2G 18300 


basic he ix-loop-helix (bHLH) family protein 


6 


AT2G23600 A 


ACL (ACETONE-CYANOHYDRIN LYASE); hydrolase 


6 


AT2G40330 


Bet v I allergen family protein 


6 


AT2G46450 


ATCNGC12 (cyclic nuceotide gated channel 12); cyclic nucleotide binding 


6 


AT3G14310* A 


ATPME3 (Arabidopsis thaliana pectin methylesterase 3) 


6 


AT3G23880* 


F-box family protein 


6 


AT3G49260* 


IQD21 (IQ- DOMAIN 21, IQ-domain 21); calmodulin binding 


6 


AT3G49670 


BAM2 (big apical meristem 2); ATP binding/protein serine/threonine kinase 




AT3G49940* 


LBD38 (LOB DOMAIN-CONTAINING PROTEIN 38) 


6 


AT4G 15390 


transferase family protein 


6 


AT4G 17870* 


similar to unknown protein [Arabidopsis tha iana] (TA R:AT4G1 7870.1) 


6 


AT4G21410* 


protein kinase fami y protein 




AT4G24780 


pectate lyase family protein 


6 


AT4G36540 


BEE2 (BR ENHANCED EXPRESSION 2); DNA binding/transcription factor 




AT4G38840 


auxin-responsive protein, putative 


6 


AT5G07580* 


DNA binding/transcription factor 


6 


AT5G14120* 


nodulin fami y protein 


6 


AT5G 14760* 


AO (L-ASPARTATE OXIDASE); L-aspartate oxidase 


6 


AT5G25460 


unknown protein 


6 


AT5G25840* 


unknown protein 


6 


AT5G39080* 


transferase family protein 


6 


AT5G61590 


AP2 domain-containing transcription factor family protein 


6 


AT5G63180 


pectate lyase family protein 


6 


AT5G64100 


peroxidase, putative 


6 


AT5G66690 


UGT72E2; UDP-glycosyltransferase/coniferyl-alcohol glucosyltransferase 


6 


AT5G66770 


scarecrow transcription factor family protein 


6 



The first column lists the AGI locus identifiers of the genes. *Genes that appear in our guard cell ABA-repressed gene set. A Genes that appear in Leonhardt et 
al.'s guard cell ABA-repressed gene set [16]. The genes that appear in our ABA-repressed gene set are marked with underlines if they are not represented on the 
Affymetrix ~8K chips used in [16]. The second column gives brief TAIR descriptions for the genes. The third column indicates the number of experiments that the 
genes are found common to. 



Wang et at. BMC Genomics 201 1, 12:216 
http://www.biomedcentral.eom/1 471 -2 1 64/1 2/216 



Page 1 2 of 24 



Table 5 Top 50 ABA-induced genes in guard cells that have not been reported in other ABA transcriptome studies 



AGI # 


Gene name 


Description 


P-value 


Fold change 


AT4G12130 




aminomethyltransferase 


1.2E-12 


6.1 


AT3G01510 




5'-AMP-activated protein kinase beta-1 subunit-related 


2.5E-12 


8.5 


AT3G15357 




unknown protein 


27E-12 


31.2 


AT5G13930 


ATCHS 


CHALGONE SYNTHASE 


47E-1 1 


261.3 


AT3G04460 


ATPEX12 


PEROXIN-12; actin binding 


63E-1 1 


3.7 


AT2G 19580 


TET2 


TETRASPANIN2 


1.1E-10 


4.3 


AT2G43018 




conserved peptide upstream open reading frame 17 


57E-10 


4.4 


AT3G14180 




transcription factor 


1 .2E-09 


3.6 


AT1 G32550 




ferredoxin family protein 


13E-09 


4.9 


AT5G 17980 




C2 domain-containing protein 


1.5E-09 


5.0 


AT5G67370 




unknown protein 


1 .8E-09 


4.5 


AT2G44660 




transferase, transferring glycosyl groups 


1 .8E-09 


3.9 


AT3G18170 




similar to unknown protein AT3G18180.1 


1 .9E-09 


3.3 


AT4G08980 




F-box family protein (FBW2) 


24E-09 


3.6 


AT1 G68470 




exostosin family protein 


24E-09 


10.3 


AT2G24150 


HHP3 


heptahelical protein 3; receptor 


2.6E-09 


1 0.0 


AT5G60790 


ATGCN1 


Arabidopsis tha iana genera contro non-repressible 1 


4.3E-09 


4.8 


AT5G61670 




heat shock protein binding/unfo ded protein binding 


5.0E-09 


3.0 


AT3G23580 


RNR2 


RIBONUCLEOTIDE REDUCTASE 2A 


6.0E-09 


6. 1 


AT4G 17420 




similar to unknown protein AT5G47420.1 


7.6E-09 


3.8 


AT1 G67960 




similar to unnamed protein product GB:CA042391.1 


8.2E-09 


2.2 


AT5G45920 




carboxylesterase 


1 .3E-08 


5.9 


AT5G 16990 




NADP-dependent oxidoreductase, putative 


1.5E-08 


2.8 


AT2G35700 


ERF38 


ERF FAMILY PROTEIN 38; transcription factor 


1 .6E-08 


2.5 


AT4G1 1370 


RHA1A 


RING-H2 finger A1A; zinc ion binding 


1 7E-08 


5.9 


AT2G04400 




indole-3-glycero phosphate synthase (IGPS) 


1 .8E-08 


3.7 


AT5G07250 


ATRBL3 


ARABIDOPSIS RHOMBOID-LIKE PROTEIN 3 


2.3E-08 


3.4 


AT3G 13040 




myb fami y transcription factor 


2.4E-08 


2.8 


AT4G00370 


ANTR2 


anion transporter 2 


2.5E-08 


2.7 


AT4G34930 




1 -phosphatidylinositol phosphodiesterase-re ated 


2.7E-08 


3.4 


AT1G01570 




fringe-related protein 


2.8E-08 


2.7 


AT1 G23750 




DNA-binding protein-related 


2.9E-09 


3.1 


AT5G63330 




DNA-binding bromodomain-containing protein 


34E-08 


3.1 


AT1G 18360 




hydrolase, alpha/beta fold family protein 


37E-08 


6.2 


AT5G55090 


MAPKKK15 


Mitogen-activated protein kinase kinase kinase 15 


4.8E-08 


5.9 


AT1 G6371 0 


CYP86A7 


cytochrome P450, fami y 86, subfamily A, polypeptide 7 


5.9E-08 


10.6 


AT2G29730 


UGT71D1 


UDP-GLUCOSYL TRANSFERASE 71 D1 


7.5E-08 


8.1 


AT5G03160 


ATP58IPK 


ARABIDOPSIS HOMOLOG OF MAMALLIAN P58IPK 


8.8E-08 


2.6 


AT3G50850 




similar to unknown protein AT5G49560.1 


1 .0E-07 


3.7 


AT1 G03630 


POR G 


PROTOCHLOROPHYLLIDE OXIDOREDUCTASE 


1.1E-07 


2.8 


AT2G39020 




GCN5-related N-acetyltransferase (GNAT) family 


1.2E-07 


3.6 


AT1G10960 


ATFD1 


FERREDOXIN 1; 2 iron, 2 sulfur cluster binding 


1.2E-07 


2.5 


AT1G19190 




hydrolase 


1 .8E-07 


6.7 


AT2G45300 




3-phosphoshikimate 1-carboxyvinyltransferase 


1 .8E-07 


2.9 


AT2G44060 




late embryogenesis abundant family protein 


1 .9E-07 


3.3 


AT4G 10730 




kinase 


2.1E-07 


2.9 


AT2G47670 




invertase/pectin methylesterase inhibitor family protein 


2.1E-07 


3.1 


AT3G54960 


ATPDIL1-3 


PDI-LIKE 1-3; thiol-disulfide exchange intermediate 


2.3E-07 


3.6 


AT4G22320 




similar to unknown protein TAIR:AT5G55210.1 


2.6E-07 


1.7 


AT5G61370 




pentatricopeptide (PPR) repeat-containing protein 


2.6E-07 


5.2 



The fold changes are calculated by averaging three biological replicates and comparing the average expression values with ABA treatment to those without ABA 
treatment. The genes with underlines also appear in our leaf ABA-induced genes. 



Wang et at. BMC Genomics 201 1, 12:216 
http://www.biomedcentral.eom/1 471 -2 1 64/1 2/216 



Page 1 3 of 24 



Table 6 Top 50 ABA-repressed genes in guard cells that have not been reported in other ABA transcriptome studies 



AG I # 


Gene name 


Description 


P-value 


Fold change 


AT2G45 1 20 




zinc finger (C2H2 type) family protein 


4.7E-15 


14.7 


AT1 G60630 




leucine-rich repeat family protein 


8.2E-13 


6.7 


AT3G 19850 




phototropic-responsive NPH3 family protein 


3.3E-12 


13.3 


ATI G1 9620 




unknown protein 


1.1E-11 


9.0 


AT1G1 1340 




S-locus lectin protein kinase family protein 


1.5E-1 1 


23.6 


AT4G04955 


ATALN 


ARABIDOPSIS ALLANTOINASE 


2.8E-1 1 


3.3 


AT3G19120 




unknown protein 


34E-1 1 


10.6 


AT4G37870 


PCK1 


PHOSPHOENOLPYRUVATE CARBOXYKINASE 1 


37E-1 1 


-1.5 


AT4G01770 


RGXT1 


RHAMNOGALACTURONAN XYLOSYLTRANSFERASE 1 


1.2E-10 


7.3 


AT5G47780 


GAUT4 


Ga acturonosyltransferase 4 


24E-1 0 


■'1.2 


AT5G52120 


ATPP2-A14 


Ph oem protein 2-A14; carbohydrate binding 


2.9E-10 


6.1 


AT1G21540 




AMP-binding protein, putative 


4.6E-10 


6.8 


AT1 G0901 0 




glycoside hydrolase family 2 protein 


7.3E-10 


3.4 


AT1G28010 


PGP14 


P-GLYGOPROTEIN 14; ATPase 


1 .2E-09 


12.6 


AT3G52870 




ca modulin-binding famiy protein 


17E-09 


10.2 


AT1 G75880 




fami y II extracellular lipase 1 (EXL1) 


2.5E-09 


-1.2 


AT2G25780 




unknown protein 


3.8E-09 


48.1 


AT2G45810 




DEAD/DEAH box helicase, putative 


6.7E-09 


-1.1 


AT3G62780 




C2 domain-containing protein 


7.6E-09 


6.7 


AT5G57580 




ca modulin-binding protein 


1 .5E-08 


2.8 


AT3G09770 




zinc finger (G3HC4-type RING finger) family protein 


1 .6E-08 


3.3 


AT4G08920 


CRY1 


GRYPTOCHROME 1 


17E-08 


■'1.8 


AT2G34160 




nucleic acid binding 


17E-08 


4.9 


AT3G24140 


FMA 


FMA (FAMA); DNA binding/transcription activator 


1 .9E-08 


5.3 


AT1 G75450 


CKX5 


CYTOKININ OXIDASE 5; cytokinin dehydrogenase 


2.1 E-08 


2.9 


AT3G1 1540 


SPY 


SPINDLY; transferase, transferring glycosy groups 


2.1 E-08 


2.6 


AT3G51390 




zinc finger (DHHG type) fami y protein 


2.3E-08 


2.8 


AT1G12730 




ce I division cycle protein-related 


2.5E-08 


5.3 


AT1G1 1260 


STP1 


SUGAR TRANSPORTER 1 


2.6E-08 


3.3 


AT2G 16640 


ATTOC132 


MULTIMERIC TRANSLOCON COMPLEX IN THE OUTER ENVELOPE MEMBRANE 132 


2.9E-08 


2.7 


AT5G53890 




leucine-rich repeat transmembrane protein kinase 


3.3E-08 


3.3 


AT5G1 1580 




UVB-resistance protein-related 


3.6E-08 


2.7 


AT2G22300 


SR1 


SIGNAL RESPONSIVE 1; calmodulin binding 


4.5 E-08 


6.5 


AT5G23880 


CPSF100 


CLEAVAGE AND POLYADENYLATION SPECIFICITY FACTOR 


4.6E-08 


1.8 


AT2G44065 




ribosomal protein L2 family protein 


5.1 E-08 


2.9 


AT3G 17020 




universal stress protein (USP) fami y protein 


5. 2 E-08 


2.7 


AT1 G74040 


IMS1 


IMS1; 2-isopropylma ate synthase 


5.4E-08 


3.7 


AT1G17230 




protein binding/protein kinase 


5. 5 E-08 


3.4 


AT3G05160 




sugar transporter, putative 


5.9E-08 


-1.8 


AT2G37420 




kinesin motor protein-re ated 


6.0E-08 


1 1.0 


AT1G63120 


ATRBL2 


ARABIDOPSIS THALIANA RHOMBOID-LIKE 2 


6.3 E-08 


3.1 


AT5G64560 




magnesium transporter CorA-like family protein 


8.8E-08 


2.4 


AT5G63620 




zinc-binding dehydrogenase family protein 


9.3E-08 


2.6 


AT4G26830 




hydrolase, hydrolyzing O-glycosyl compounds 


9.9E-08 


7.3 


ATI G2 1400 




2-oxoisovalerate dehydrogenase, putative 


1.3E-07 


3.0 


AT5G48950 




thioesterase family protein 


1.3E-07 


5.2 


AT5G37450 




leucine-rich repeat transmembrane protein kinase 


1 4E-07 


9.9 


AT3G05660 


AtRLP33 


Receptor Like Protein 33; kinase/protein binding 


1.5E-07 


3.1 


AT3G58840 




similar to myosin heavy chain-related AT1G06530.1 


1 .8E-07 


3.9 


AT4G38300 




glycosyl hydrolase family 10 protein 


2.1E-07 


3.8 



The fold changes are calculated by averaging three biological replicates and comparing the average expression values with ABA treatment to those without ABA 
treatment. The genes with underlines also appear in our leaf ABA-repressed genes. 
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Response to abiotic stimulus (92, 79) 
Response to endogenous simulus {67, 55) 
Response to stress (106, 107) 
Lipid metabolic process (36, 26) 
Catalytic activity (241 , 200) 
Secondary metabolic process (20, 23) 
Vacuole (28, N/A) 
Transcription factor activity (66, N/A) 
Lipid binding (N/A, 12) 
Translation (3, 0) 
RNA binding (5, 0) 
Nucleobase metabolic process (11,3) 
Protein metabolic process (N/A, 24) 
Nucleotide binding (N/A, 22) 
Structural molecule activity (N/A, 1) 




I GC 
I LF 



underrepresentation overrepresention 
-log 10 (P-value) 



B 

Plastid (81 , 85) 
Transferase activity (80, 79) 
Plasma membrane (53, 59) 
Nucleotide binding (60, 75) ■ 
Response to abiotic stimulus (51 , N/A) ■ 
Response to stress (54, N/A) 
Transporter activity (36, N/A) 
Endoplasmic reticulum (16, N/A) 
Biosynthetic process (67, N/A) 
Protein modification process (N/A, 39) ■ 
Nucleobase metabolic process (N/A, 35) ■ 
Translation (6, N/A) 



underrepresentation overrepresention 
-log 1D (P-value) 

Figure 7 Significantly enriched functional categories in guard 
cell and leaf ABA-responsive genes (A) ABA-induced genes. (B) 
ABA-repressed genes. Numbers in parentheses indicate the numbers 
of genes annotated with the corresponding functional categories in 
guard cell and leaf ABA-regulated gene sets respectively. 'N/A' 
means the enrichment of the function is not significant. 



Cross-regulation by ABA and heterotrimeric G proteins 

Heterotrimeric G-proteins, composed of a, P and y sub- 
units, participate in a wide range of crucial signaling 
pathways in eukaryotes [58,59]. ABA signaling is known 
to interact with heterotrimeric G-protein signaling in 
both developmental and stress responses [60-63]. In a 
previous study, we found that ABA signaling also cross- 
talks with G protein signaling at the level of the tran- 
scriptome [31]. In the present study, we used the gen- 
ome-wide guard cell microarray data from wild-type, 
Ga (gpal), G|3 (agbl) and Gy double {agbl gpal) 
mutant plants, with and without ABA treatment to sys- 
tematically investigate ABA regulation of gene expres- 
sion in guard cells and compared it with the ABA- 
regulated transcriptome of leaves from the same geno- 
types. We identified a number of G-protein-independent 
ABA-regulated genes and G-protein-dependent ABA- 
regulated genes, which will facilitate the screening and 
identification of novel guard cell signaling genes. Figure 
9 demonstrates the heat map of expression patterns of 
these genes (drawn by Matrix2png [64]) and shows that 
in both guard cells and leaves the number of G-protein- 
independent ABA-regulated genes is dominant, consis- 
tent with the fact that ABA mediates many stress and 
developmental signaling pathways without any known 
participation from the heterotrimeric G protein. The 
number of G-protein-dependent ABA-regulated genes in 
leaves is much larger than in guard cells, suggesting that 
the G-protein may be more involved in ABA signaling 
in mesophyll cells than in guard cells. Six (out of 9) 
guard cell G-protein-dependent ABA-induced genes and 
14 (out of 58) leaf G-protein-dependent ABA-induced 
genes have ABREs in their promoter regions, which is 
statistically significant. 



Guard cell signaling and development genes regulated at 
the transcript level by ABA 

A number of genes and proteins in guard cells have 
been revealed to function in stomatal movements in 
response to stimuli, including ABA, CO2, light and 
pathogens. In [55], a total of 67 guard cell signaling 
genes was compiled from published literature. Li and 
Assmann [56] reviewed selected guard cell signaling and 
development genes. In addition, Kwak et al. collected 69 
guard cell signaling genes [13]. We combined these sets 
as well as 22 additional genes with functional guard cell 
roles as described by Kim et al. [14] and 24 genes 
involved in guard cell patterning and development as 
described by Dong and Bergmann [57], to obtain a list 
of 149 known guard cell signaling or development 
genes/proteins (Additional file 6). Of these guard cell 
components, 29 are significantly regulated by ABA at 
the transcript level (Figure 8). 



Cross-regulation by ABA and other hormones 

The phenomenon of cross-talk between plant hormones 
is well-established at the physiological level [10,65,66]. 
Nemhauser et al. [22] performed an analysis of seven 
hormone-treated transcriptomes generated by Goda et 
al. [67] from 7-day-old Arabidopsis seedlings and identi- 
fied target genes regulated by the hormones abscisic 
acid (ABA), gibberellic acid 3 (GA), indole-3-acetic acid 
(IAA, auxin), 1-amino-cyclopropane-l-carboxylic acid 
(ACC, ethylene precursor), zeatin (CK, cytokinin), bras- 
sinolide (BL, brassinosteroid), and methyl jasmonate 
(MJ, jasmonate). Six of these seven hormones both posi- 
tively and negatively regulate significant numbers of 
genes at the transcript level, identified by stringent ana- 
lysis using two different methods [22], whereas the 
seventh hormone, GA, regulates only a few genes identi- 
fied by a low- stringency linear model [22]. We took the 
intersection of the hormone-regulated gene sets identi- 
fied by the rank product and linear model methods used 



Wang ef al. BMC Genomics 201 1, 12:216 
http://www.biomedcentral.eom/1 471 -2 1 64/1 2/216 



Page 1 5 of 24 



-4 



T 



-2 0 2 

log 2 (fold change) 



AT5G67030 ABA1 (ABA DEFICIENT 1 ); zeaxanthin epoxidase 

AT4G34000 ABF3 (ABSCISIC ACID RESPONSIVE ELEMENTS-BINDING FACTOR 3) 

AT3G19290 ABF4 (ABRE BINDING FACTOR 4); DNA binding 

AT4G26080 ABM (ABA INSENSITIVE 1); calcium ion binding; protein serine/threomine phosphatase 

AT5G57050 ABI2 (ABA INSENSITIVE 2); protein serine/threonine phosphatase 

AT5G67300 ATMYB44/ATMYBR1/MYBR1 (MYB DOMAIN PROTEIN 44); DNA binding 

AT3G11410 ATPP2CA (ARABIDOPSIS THALIANA PROTEIN PHOSPHATASE 2CA) 

AT1G72770 HAB1 (HOMOLOGY TO ABI1); protein serine/threonine phosphatase 

AT2G46370 JAR1 (JASMONATE RESISTANT 1) 

AT1G01360 PYL9/RCAR1 (regulatory component of ABA receptor 1 ) 

AT1G28010 ABCB14/PGP14 (P-GLYCOPROTEIN 14); ATPase 

AT3G53720 ATCHX20 (CATION/H* EXCHANGER 20); monovalent cation:proton antipor 

AT5G40890 ATCLC-A (CHLORIDE CHANNEL A); anion channel/voltage-gated chlorid (tonoplast) 

AT3G45640 ATMPK3 (MITOGEN-ACTIVATED PROTEIN KINASE 3) 

AT1G12110 CHL1/NRT1 .1 (NITRATE TRANSPORTER 1 .1 ); transporter 

AT4G08920 CRY1 (CRYPTOCHROME 1) 

AT3G24140 FMA (FAMA); DNA binding/transcription activator 

AT5G46240 KAT1 (K* ATPASE 1); inward rectifier potassium channel 

AT4G18290 KAT2 (K + ATPase 2); inward rectifier potassium channel 

AT1G08810 MYB60 (myb domain protein 60); DNA binding/transcription factor 

AT1G37130 NIA2 (NITRATE REDUCTASE 2) 

AT2G18960 OST2/AHA1 (ARABIDOPSIS H + ATPASE 1); ATPase 

AT5G46790 PYL1 (PYR1-like 1); regulatory components of ABA receptor 12 

AT2G26040 PYL2 (PYR1-like 2); regulatory components of ABA receptor 14 

AT2G38310 PYL4 (PYR1-like 4); regulatory components of ABA receptor 10 

AT5G05440 PYL5 (PYR1-like 5); regulatory component of ABA receptor 8 

AT4G17870 PYR1 (PYRABACTIN RESISTANCE 1); regulatory component of ABA receptor 

AT2G30520 RPT2 (ROOT PHOTOTROPISM 2); protein binding 

AT1G11260 STP1 (SUGAR TRANSPORTER 1); carbohydrate transmembrane transporter 



Figure 8 Known guard cell signaling and development genes regulated by ABA at the transcript level The x axis denotes the log2 ratio 
of average gene expression with ABA treatment over that without ABA treatment. The y axis lists known guard cell signaling and development 
genes. The dashed lines indicate 2-fold change. All genes depicted showed significant regulation of gene expression based on FDR-corrected P- 
value of the linear model method, and also passed the correlation threshold of the Boolean method. 



in [22] and compared our ABA-regulated genes in guard 
cells and leaves with these hormonally regulated genes. 
The overlap of our ABA-regulated genes with hormone- 
regulated gene sets was calculated both in the same 
direction, that is, genes up- (down-) regulated by a hor- 
mone are compared with genes up- (down-) regulated 
by another hormone, and in the antagonistic direction. 
Representation factor and P-values calculated by the 
hypergeometric distribution (see Methods) were used to 
evaluate the overlaps. 

Figure 10 shows the significance of the overlaps of our 
ABA-regulated genes in guard cells and leaves with the 
six hormone-regulated gene sets of Nemhauser et al. 
[22]. Our ABA-regulated genes in both guard cells and 
leaves have highly significant overlap with ABA-regu- 
lated and MJ-regulated genes in [22]. ABA-induced 



genes have a larger overlap with MJ-induced genes in 
leaves than in guard cells. Although the numbers are 
smaller, guard cell ABA-induced genes also have a sig- 
nificant overlap with IAA-induced genes, and leaf ABA- 
induced genes have a significant overlap with ACC- 
induced genes. Guard cell ABA-repressed genes have a 
significant overlap with BL-repressed genes. Overlap of 
our ABA-regulated genes with hormone-regulated gene 
sets in an antagonistic way, although not significant, is 
observed in each hormone-regulated gene set (Table S4 
in Additional file 3). In some instances, the extent of 
antagonistic overlap is similar to or even larger than 
that corresponding to the comparison in the same direc- 
tion. For example, ABA-repressed genes in guard cells 
and leaves have 10 and 8 genes in common with IAA- 
induced genes, and ABA-induced genes in guard cells 
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d b gpalagbl wt d b gpalagbl wt 



d b gpalagbl wt d b gpalagbl wt 




G-protein-independent 
ABA-regulated genes 




G-protein-independent 
ABA-regulated genes 



G-protein-dependent 
ABA-regulated genes 



J G-protein-dependent 
ABA-regulated genes 

Figure 9 Expression patterns of ABA-responsive genes in guard cells and leaves (A) Expression patterns of ABA-regulated genes in guard 
cells. (B) Expression patterns of ABA-regulated genes in leaves. The rows of each pattern correspond to genes, and the columns correspond to 
conditions, i.e. four genotypes, agbl gpal double mutant (db), gpal mutant, agb\ mutant, wild-type (wt), without or with ABA treatment. 



and leaves have 14 and 13 common genes with MJ- 
repressed genes. Nemhauser et al. defined marker genes 
as those genes specifically regulated by one hormone in 
a high stringency analysis and not by any other hor- 
mone even in a lower stringency analysis [22]. However, 
some genes designated as marker genes specific for hor- 
mones other than ABA are actually found to be ABA- 



regulated in our guard cell or leaf transcriptomes (Table 
S5 and Table S6 in Additional file 3). 

Discussion 

Our identified ABA-regulated guard cell and leaf tran- 
scriptomes are supported not only by Q-PCR verifica- 
tion but also by their significant overlap with the ABA- 




ABA ACC BL IAA MJ CK 




ABA ACC BL IAA MJ CK 



B 



Figure 10 Comparison of our ABA-responsive genes with the hormone-regulated genes in [22] GC: guard cell; LF: leaf. (A) Hormone- 
induced genes. The numbers of hormone-induced genes in [22] are: ABA 512, ACC 34, BL 28, IAA 198, MJ 522, and CK 60. (B) Hormone- 
repressed genes. The numbers of hormone-repressed genes in [22] are: ABA 270, ACC 23, BL 33, IAA 61, MJ 231, and CK 16. The dashed lines 
indicate the threshold for significance (P < 0.01). The numbers over the bars represent the numbers of our ABA-regulated genes that overlap 
with the hormone-regulated genes in [22]. 
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regulated transcriptomes of other tissues with respect to 
the promoter motifs that are utilized, the genes that are 
regulated, and the functional roles of the encoded gene 
products in stress tolerance as ascertained from GO 
analysis. However, as discussed in detail below, a deeper 
study of these new transcriptomes reveals many interest- 
ing features. 

Promoters 

As expected, known promoter motifs conferring ABA- 
and stress-related gene regulation, such as ABREs and 
LTREs, are over-represented in the promoters of both 
our guard cell and leaf ABA-induced gene sets. How- 
ever, promoters of ABA-regulated genes of guard cells 
but not leaves are enriched in MYB and DRE/CRT ele- 
ments. These differences suggest differential utilization 
of specific transcription factors in guard cells vs. leaves. 
Functional support for this hypothesis is offered by pub- 
lished analyses of MYB family members and stomatal 
regulation. Two members of the Arabidopsis R2R3-MYB 
family, MYB 15 and AtMYB44, are expressed in guard 
cells and their overexpression lines show enhancement 
of ABA-induced stomatal closure and drought tolerance 
[68,69]. Comparable functional data detailing guard cell 
roles for the DRE/CRT binding factors are presently 
lacking; our analysis of promoter motifs suggests that 
DRE binding proteins will also be found to participate 
in guard cell ABA responses. Analyses in other tissues 
have shown that genes with DRE elements in their pro- 
moters are regulated by dehydration but not by ABA, 
although ABA-sensitivity may be conferred by the con- 
certed action of the DREB and ABF proteins [70]. When 
our guard cell ABA-regulated gene set is assessed for its 
complement of known transcription factors, we find that 
DREB1A, ABF3 and ABF4 show concerted ABA-induc- 
tion in guard cells. 

We hypothesized that, given the specialized nature of 
guard cell responses to ABA, these cells might also uti- 
lize specialized promoter motifs. One motif, GTCGG, 
which is not a subsequence of any known motifs, was 
thereby provisionally identified, awaiting functional veri- 
fication. In addition, a greater number of different 
enriched motifs are found for guard cell ABA-induced 
genes than for leaf ABA-induced genes, suggesting that 
transcription factors in guard cells may target a greater 
diversity of cw-acting regulatory elements in ABA 
responses as compared to those in leaves. 

Only a few motifs, including CAAGTTG, emerged as 
significant from the promoters of ABA-repressed genes. 
For both guard cells and leaves, there are more signifi- 
cant 5-10mer motifs in ABA-induced genes than in 
ABA-repressed genes. Since the numbers of ABA- 
induced and ABA-repressed genes are of similar magni- 
tude in our datasets, this result suggests that evolution 



may have favored a greater diversity in mechanisms for 
regulation of gene repression than for gene induction, 
thus is it more difficult for any one ABA-repression 
motif to achieve statistical significance. 

Comparison with other ABA transcriptome studies 

The fact that many classic ABA- and stress-related pro- 
moter motifs are enriched upstream of genes in our 
ABA-regulated gene sets led us to predict that many of 
the genes identified by our analyses would also be pre- 
sent in comparable transcriptomes derived from other 
tissues. Indeed, we find over half of all ABA-regulated 
genes in our guard cell and leaf datasets were also iden- 
tified in other studies, with the highest percentage over- 
lap (86.6%) found for ABA-induced genes in our leaf 
gene set, probably reflecting the fact that leaf tissue is 
included as a component tissue (e.g. as a component of 
seedlings and whole plants) in several previous studies. 
On the other hand, it is interesting that there is only a 
single ABA-regulated gene, AT1G01470: LEA14 (LATE 
EMBRYOGENESIS ABUNDANT 14) found in common 
to all ABA transcriptome studies that we analyzed. This 
result hints that there might be greater flexibility in how 
tissues achieve a common end result of stress tolerance 
than in how they perform other functions such as pri- 
mary metabolism. 

Despite the lack of genes found in common to all stu- 
dies, we nevertheless identified core ABA-regulated gene 
sets. Many of these core ABA-induced genes have 
already been assigned functional roles in response to 
ABA or another stress, or are members of a family (e.g. 
PP2C phosphatases [33-35]) already known to play such 
a role. In contrast, fewer ABA-repressed genes (Table 4) 
have known relationships to stress signaling or toler- 
ance. Taken together with the observation that it has 
been more difficult to identify promoter motifs for 
ABA-repression than for ABA-induction, and with the 
fact that we had to reduce the number of experiments 
to 6 to identify a shared set of -50 ABA-repressed 
genes, these results suggest a greater between-tissue 
diversity in the roles of ABA-repressed genes than those 
of ABA-induced genes. For example, our guard cell 
dataset includes 56 of the 67 core ABA-induced genes 
but only 24 of the 47 core ABA-repressed genes. 

During drought stress, it is optimal for guard cells to 
lose water to promote stomatal closure while it is opti- 
mal for other cell types to retain water to maintain cel- 
lular hydration. We therefore expected to find some 
ABA-regulated genes unique to guard cells, and this 
expectation was satisfied. We identified -150 ABA- 
induced and -150 ABA-repressed genes in guard cells 
that were not identified as ABA-regulated genes in any 
of 11 previous studies of ABA-regulated transcriptomes 
(see Table S3 in Additional file 3). These genes 



Wang et al. BMC Genomics 201 1, 12:216 
http://www.biomedcentral.eom/1 471 -2 1 64/1 2/216 



Page 1 8 of 24 



represent potential guard cell-specific ABA-regulated 
genes. Consistent with this hypothesis, our predicted 
guard cell-specific motif GTCGG is significantly 
enriched in the upstream regions of these genes (P < 
1.2E-04). Supporting the idea that these genes may have 
novel, as yet undescribed ABA-related functions in 
guard cells, many fewer genes in the lists of new and 
most significantly ABA-regulated guard cell genes 
(Tables 5, 6) are annotated as ABA- or stress-related 
than in the lists of core ABA-regulated genes (Tables 3, 
4). 

Additional file 4 provides the full list of the newly 
identified ABA-responsive genes of guard cells. Some of 
these genes already have been reported to be involved in 
stomatal movements or guard cell signaling. JAR1 
(AT2G46370) is Jasmonate Resistant 1 and has been 
detected in guard cells [16]. This gene product has been 
reported to enhance ROS production and ABA sensitiv- 
ity of guard cells [71] and its expression is ABA-induced 
in our guard cell data. Other hormone-regulated genes, 
such as the IAA inducible gene IAA16 (AT3G04730) 
and cytokinin response factor 3 CRF3 (AT5G53290), are 
also regulated by ABA in guard cells according to our 
microarray analysis. In addition, several Ca 2+ -related 
genes are identified as new ABA-induced genes in guard 
cells, such as AtCAMBP25 (AT2G41010) and CBL2 
(AT5G55990). CRY1/BLU1 (AT4G08920) encodes the 
blue light photoreceptor cryptochrome 1 and has been 
reported to participate in blue light-regulation of stoma- 
tal movements [72], while NPH3 functions downstream 
of the PHOT1/PHOT2 blue light photoreceptors which 
mediate blue light-specific stomatal opening [73]. Our 
results indicate that these genes are ABA-repressed in 
guard cells, which might contribute to the phenomenon 
of ABA-repression of light-induced stomatal opening, 
particularly over a longer time scale. Sucrose is a signifi- 
cant osmoticum maintaining stomatal opening in the 
afternoon [74], and the sugar transporter STP1 has 
already been described as a guard cell-specific transpor- 
ter that shows a peak in expression in the afternoon 
[75], in concert with the daily timing of sucrose accu- 
mulation into guard cells [74]. The observed repression 
of STP1 and the putative sugar transporter AT3G05160 
by ABA might also contribute to long-term inhibition of 
stomatal opening. 

GO functional analysis 

The GO category lipid metabolic processes is over-repre- 
sented in ABA-induced genes of guard cells and leaves. 
Little is known concerning the relationship between 
ABA and lipid signaling in mesophyll cells; however, the 
guard cell result is consistent with functional evidence 
for roles of inositol phosphates (produced by the action 
of phospholipase C) and phosphatidic acid (produced by 



the action of phospholipase D) in ABA inhibition of sto- 
matal opening [76-79]. In fact, among ABA-induced 
genes in this category we find AT5G58700, a phosphoi- 
nositide-specific phospholipase C family protein, and 
AT2G22240, an inositol-3-phosphate synthase, as well 
as two PLD isoforms: PLDS (AT4G35790; also con- 
firmed in our Q-PCR analysis) and PLDC (AT3G05630), 
suggesting that activation of these genes may contribute 
toward long-term suppression of stomatal opening 
under drought conditions. 

Conversely, several functional categories are underre- 
presented within the ABA-induced gene sets, including 
translation, RNA binding, protein metabolism, and 
nucleobase metabolism, consistent with the supposition 
that metabolism is generally down-regulated under 
stress conditions [80]. However, such down-regulation 
may be greater for leaves than for guard cells, given that 
GO categories associated with several types of metabolic 
processes are over-represented in leaf but not in guard 
cell ABA-repressed gene sets. This difference may reflect 
a requirement for guard cells to maintain metabolic 
activity and active stomatal regulation under drought 
conditions. One interesting observation is that the GO 
category response to abiotic stimulus is over-represented 
in both ABA-induced and ABA-repressed guard cell 
genes. Many of these genes are genes of known function 
in guard cells, and are discussed in the next section. 

Genes of known guard cell function, G protein regulation, 
and hormone cross-talk 

Our analysis reveals that many genes which encode pro- 
teins with known functions in guard cell physiology are 
also ABA-regulated at the transcript level. This phenom- 
enon was reported less robustly in a previous analysis 
[16] which queried only -1/3 of the Arabidopsis genes 
assayed here. Changes in rates of solute transport result- 
ing from non-transcription-related regulation of ion 
channels and transporters are central to rapid, osmoti- 
cally-driven control of stomatal apertures [81,82]. How- 
ever, our data show that many of the relevant transport 
proteins are regulated at the transcript level as well. 
Transcripts encoding the H + ATPase OST2, the inward 
K + channels KAT1 and KAT2, the sugar transporter 
STP1, the endosomal Na + /H + antiporter CHX20, and 
the nitrate importer CHL1 are all down-regulated in 
guard cells following ABA treatment, suggesting that 
their downregulation may be an important aspect of 
longer-duration inhibition of stomatal opening by ABA. 
Two known genes related to guard cell development, 
FAMA {FMA) which regulates proliferation of stomatal 
precursors as well as differentiation of guard mother 
cells [83], and ATMPK3, an environmentally responsive 
mitogen-activated protein kinase mediating stomatal 
development and patterning [84], are downregulated, as 
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are several signaling proteins. Several other genes, func- 
tioning in the guard cell response to C0 2 (HT2, HIC) 
and pathogens (e.g. FLS2), are slightly repressed by ABA 
(not shown), although failing to meet the threshold of 
significance used to categorize ABA-responsive genes in 
our microarray analyses. 

ABA up-regulated genes include zeaxanthin epoxidase 
(ABA1), required for ABA synthesis. In addition, up-reg- 
ulation of ABA-related transcription factors ABF3, 
ABF4, and AtMYB44 could provide positive feedback to 
the ABA response by promoting transcription of stress- 
related genes. A number of genes encoding known 
guard cell signaling proteins are up-regulated, including 
four protein phosphatase 2C's: PP2Ca, HAB1, ABI2, and 
ABI1. Since PP2C phosphatases are negative regulators 
of ABA signaling, their ABA-induced expression sug- 
gests instigation of negative feedback on the ABA 
response. This supposition is supported by the converse 
observation that five out of six genes encoding PYR/ 
PYL/RCAR type soluble ABA receptors with documen- 
ted function in guard cells [85-87] are down-regulated 
by ABA (out of the 14 receptors in PYR/PYL/RCAR 
family, only ten are represented on the ATH1 chip), as 
are two other genes encoding positive transducers of the 
ABA signal, the MAPK MPK3, and an enzyme involved 
in nitric oxide production, NIA2. 

Analysis of T-DNA insertional mutants has implicated 
heterotrimeric G proteins in ABA responses of diverse 
tissues [58,59,88]. Mutant plants with G-protein subunit 
knockouts are ABA hyposensitive in aspects of guard 
cell ion channel and stomatal regulation [60,63,89,90] 
but, unexpectedly, ABA hypersensitive in inhibition of 
seed germination and root growth [61,91], indicating 
system-specificity of G-protein effectors. Interestingly, 
the leaf transcriptome actually has more genes co-regu- 
lated by ABA and G-proteins than the guard cell tran- 
scriptome. This result suggests important but as yet 
unidentified functions of ABA/G-protein co-regulation 
in mesophyll cell physiology. Conversely, in guard cells, 
this result implies that non-ABA related functions of G- 
proteins should be investigated. Indeed, based on phe- 
notypes of gpal mutants, G-proteins also participate in 
C0 2 and pathogen responses of guard cells [92,93]. 

Our results indicate that signaling cross-talk also 
occurs at the transcript level between ABA and other 
plant hormones. ABA-regulated and MJ- regulated genes 
identified by Nemhauser et al. [22] each have highly sig- 
nificant overlap with our ABA-regulated genes in both 
guard cells and leaves, and this is consistent with ABA- 
MJ crosstalk at the physiological level, in both guard 
cells [71,94-96] and other tissues [97]. In addition, a few 
genes designated as marker genes for hormones other 
than ABA in [22] are found to be significantly ABA- 
regulated in our transcriptomes. It is possible that the 



previous analysis [22] was confounded by the use of 
whole Arabidopsis seedlings as the material for microar- 
ray analysis, which may have masked hormonal co-regu- 
lation occurring in specific tissue or cell types. Although 
our evaluation of overlap is limited to one hormone 
(ABA), these results nevertheless indicate that the con- 
clusion of [22], that plant hormones do not regulate 
common gene sets, merits re-assessment. Our work pro- 
vides a starting point to re-assess core hormone regula- 
tory modules and hormone marker genes. 

Conclusions 

This microarray analysis investigates ABA regulation of 
gene expression in guard cells in a systematic genome- 
wide manner. A number of ABA-regulated genes of 
guard cells identified here overlap with ABA-regulated 
genes of other tissues while a subset of them show 
ABA-regulation unique to this cell type. A unique cis- 
acting motif, GTCGG, associated with ABA-induction of 
gene expression specifically in guard cells, was identified 
at the bioinformatic level. Many of the genes known to 
encode ion transporters associated with stomatal open- 
ing are down-regulated by ABA, providing one mechan- 
ism for long-term maintenance of stomatal closure 
during drought. We also found examples of both nega- 
tive and positive feedback in the transcriptional regula- 
tion by ABA of known ABA-signaling genes, particularly 
with regard to the PYR/PYL/RCAR class of soluble ABA 
receptors and their downstream targets, the type 2C 
protein phosphatases. In conclusion, the results of this 
study engender new insights into the basic cell biology 
of guard cells, reveal common and unique elements of 
ABA-regulation of gene expression in guard cells, and 
set the stage for targeted biotechnological manipulations 
to improve plant water use efficiency. 

Methods 

Microarray hybridizations and Q-PCR 

Our microarray data utilized the Affymetrix ATH1 chip 
and are from hybridizations described in detail in [31]. 
All of these data are available in the Gene Expression 
Omnibus (GEO) database http://www.ncbi.nlm.nih.gov/ 
geo with accession no. GSE19520. To reprise, epidermal 
peels with guard cells as the only intact cell type were 
used as the source of RNA for our guard cell microar- 
rays [31] (Figure 1A and Figure IB), and leaves from the 
first and second rosette layer of five-week old plants 
were the RNA source for our leaf microarrays [31]. ABA 
treatments were 50 uM x 3 hr. and EtOH was used as 
the solvent control. If an epidermal peel preparation did 
not pass quality control parameters for guard cell purity 
it was not used for RNA isolation, and if an RNA pre- 
paration was not of high quality as indicated by Bioana- 
lyzer profiling, it was discarded. Four genotypes were 
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evaluated: Col (wild-type), gpal-4 mutant, agbl-2 
mutant, and agbl-2 gpal-4 double mutant, where GPA1 
(AT2G26300) and AGB1 (AT4G34460) are the Arabi- 
dopsis G protein a subunit and |3 subunit respectively. 
For each type of sample (guard cells or leaves), three 
independent biological replicates were obtained, result- 
ing in a total of 48 microarray hybridizations (2 sample 
types x 4 genotypes x two treatments x 3 replicates). 

For Q-PCR analysis, cDNA was aliquoted and kept at 
4°C throughout each Q-PCR experiment to avoid discre- 
pancy in the data because of freeze-thaw cycles. Real 
time PCR was performed using pre-mix containing 
SYBR-Green intercalating dye (BioRad). Actin was used 
as an internal control {Actin 2: AT3G 18780 and Actin 8: 
AT1G49240) [98]. The positions of the oligonucleotide 
primers used for real-time PCR was chosen so that the 
size of all PCR products was between 100 and 150-bp. 
The suitability of the oligonucleotide sequences in terms 
of efficiency of annealing was evaluated in advance 
using the Primer 3 program. Q-PCR experiments were 
repeated thrice independently, and the data were aver- 
aged. The data obtained were analyzed with IQ5 soft- 
ware (Bio-Rad). 

Identification of differentially expressed genes 

Our ABA-regulated gene sets consisted of the intersect- 
ing set of genes identified by two methods: a Boolean 
method and a linear model method (Figure 11). 

In [31], we adopted a Boolean framework F(ABA, 
GPA1.AGB1) = B(ABA, A(GPA1, AGB1)) + C AB a to 
describe the regulatory modes of ABA and the G 



protein and identify ABA- and/or G protein-regulated 
genes. In this framework, GPA1 and AGB1 are Boolean 
variables which can have two states: 1 denoting 'on' (not 
knocked out) and 0 denoting 'off (knocked out). ABA is 
also a Boolean variable: ABA = 1 indicates the presence 
of ABA (i.e. ABA treatment) and ABA = 0 represents 
the absence of ABA (i.e. solvent control). According to 
this framework, B 4 (ABA, A) = ABA and B 13 (ABA,A) = 
not ABA determine two Boolean rules describing ABA- 
regulation of gene expression independent of the hetero- 
trimeric G protein. The rest of the Boolean rules B 
(ABA, A) describe G-protein-only regulation or the co- 
regulation of ABA and the G protein on gene expres- 
sion. A correlation measure is used to assign differen- 
tially expressed genes to these Boolean rules [31]. In the 
present study, we used 1.5 as the correlation score 
threshold for ABA-only regulated genes to designate 
them as G-protein-independent ABA-regulated genes. 
This threshold controls the false discovery rate (FDR) 
within 0.001. In addition, those ABA-G-protein co-regu- 
lated genes that also exhibit greater than 2-fold expres- 
sion change (with FDR less than 0.05) in response to 
ABA in the wild-type are designated as G-protein- 
dependent ABA regulated genes in the present study. In 
other words, all ABA-G-protein co-regulated genes 
described here satisfy the condition of showing ABA- 
regulation in the wild-type Col genotype. 

To make the identification of differentially expressed 
genes with respect to ABA more reliable, we also used 
linear models combined with empirical Bayes methods 
to determine ABA-responsive genes [32]. This method 
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Figure 11 Overlap distribution of ABA-responsive genes obtained by the two different methods. (A) Guard cells; (B) Leaves "BM: ABA" 
and "Limma: ABA" denote the G-protein-independent ABA-regulated genes obtained by the Boolean method and the linear model method, 
respectively. "BM: ABA-Gp" and "Limma: ABA-Gp" denote the G-protein-dependent ABA-regulated genes obtained by the Boolean method and 
the linear model method, respectively. The numbers with underlines are the total numbers of genes in the corresponding categories. For guard 
cells, the ABA-regulated gene set utilized here is the sum of the 1059 + 90 + 7 + 17 genes. For leaves, the ABA-regulated gene set utilized here 
is the sum of the 637 + 1 90 + 4 + 206 genes. 
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is implemented by the limma package embedded in the 
Bioconductor project [99]. We followed the 2x4 fac- 
torial design process, using +ABA/-ABA as one factor 
and genotypes agbl gpal mutant, gpal mutant, agbl 
mutant, and wild-type as another factor. Contrasts of 
interest (ABA treated wild-type versus control wild-type, 
control mutants versus control wild-type, ABA treated 
mutants versus ABA treated wild-type) were extracted 
and the P-values for moderate f-tests were adjusted by 
Benjamini and Hochberg's method to control for the 
false discovery rate (FDR). We chose FDR <0.001 as a 
cut-off to select significantly differentially expressed 
genes. The P-value threshold for the contrast of ABA 
treated wild-type versus control wild-type was set at 
0.0001, and the P-value threshold for the contrast of 
genotypes versus wild-type was set to 0.01. Those genes 
that are only significantly differentially expressed with 
respect to ABA are designated G-protein-independent 
ABA-regulated genes. Those genes that are significantly 
differentially expressed with respect to both ABA and 
genotypes are called G-protein-dependent ABA-regu- 
lated genes. As with the Boolean method, all genes 
designated as ABA-G-protein co-regulated from the lin- 
ear model approach satisfy the condition of showing 
ABA-regulation in the wild-type Col genotype. 

The Boolean method identified 1323 G-protein-inde- 
pendent ABA-regulated genes and 32 G-protein-depen- 
dent ABA-regulated genes in guard cells, with 1024 and 
408 genes, respectively, identified in leaves. The linear 
model method identified 1379 G-protein-independent 
ABA regulated genes and 204 G-protein-dependent 
ABA-regulated genes in guard cells, with 1040 and 998 
genes, determined in leaves. The G-protein-independent 
ABA-regulated genes determined by the two methods 
(intersection of green and red circles in Figure 11) and 
the ABA-regulated genes identified by both methods but 
with inconsistency in the existence of G-protein regula- 
tion (intersection of green and yellow circles plus inter- 
section of red and blue circles in Figure 11) are 
designated G-protein-independent ABA-regulated genes 
in our final gene list. G-protein-dependent ABA-regu- 
lated genes that are confirmed by both methods (inter- 
section of yellow and blue circles in Figure 11) are 
reported as G-protein-dependent ABA-regulated genes 
in our final gene list. The genes in each subsector of 
Figure 11A and Figure 11B are given in Additional file 7 
and Additional file 8, respectively. The genes that are 
present in guard cells or leaves but not regulated by 
ABA at the transcript level (i.e. not present in any of 
the 8 subsectors of Figure 11A and Figure 11B) are 
given in Additional file 9 and Additional file 10, 
respectively. 



Promoter motif analysis 

We performed a statistical analysis of the promoter 
sequences of ABA-responsive genes in guard cells and 
leaves. The 1000-bp upstream regions of all Arabidopsis 
gene sequences were obtained from the FTP site of the 
Arabidopsis Information Resource (TAIR) (TAIR9_up- 
stream_1000_20090619.txt) [54]. Motif search was per- 
formed by writing custom Python scripts http://www. 
python.org which search a given motif pattern with an 
exact match without insert or mismatch. The signifi- 
cance of motif enrichment was computed by the hyper- 
geometric cumulative distribution function. Specifically, 
when we examined the enrichment of a known czs-regu- 
latory element, we counted the number of ABA-regu- 
lated genes in guard cells or leaves and the number of 
all genes in the chips whose promoter regions contain 
the cw-regulatory element. A P-value was obtained 
according to occurrences of this cw-regulatory element 
in ABA-regulated genes and in all genes. Positional dis- 
tribution of a cw-regulatory element in ABA-regulated 
genes was obtained by shifting a window of the same 
length to scan the promoter region of each gene, and 
recording the occurrence positions. To identify possible 
new motifs that may be involved in ABA-regulated gene 
expression, we generated sets of 5-10mer sequences and 
calculated the enrichment significance of each sequence 
in our ABA-regulated gene sets. As there are many 
more significant motifs in ABA-induced genes than in 
ABA-repressed genes, we set 1.0E-10 and 1.0E-04 as the 
significance thresholds for new motifs of ABA-induction 
and ABA-repression respectively. 

Comparison of hormone-regulated gene sets 

In the comparison of our ABA-regulated genes with 
previous ABA/hormone transcriptome studies, we used 
a representation factor (RF) and its associated P-value to 
evaluate whether the overlap of two gene sets from the 
same background is significant or not. The representa- 
tion factor is defined by the ratio of the size of the real 
overlap to the expected number of common genes 
between the two gene sets: 



RF(A, B) 



\a n B| 

N 

N X N 



\A n B\N 
\A\\B\ ' 



where \A\ is the number of genes in the gene set A, \ 
B\ is the number of genes in the gene set B, \A(~\B\ 
represents the number of genes common to gene set A 
and gene set B, and N is the total number of genes in 
the Arabidopsis genome (in this study, it is the number 
of genes covered by the Affymetrix ATH1 gene chip). 
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The P-value associated with the representation factor 
was calculated by the hypergeometric distribution: 



P(X = k) = 



f\A\\ (N-\A\\ 
\ k ) \ \B\-k) 



where k = |AnB|, and the brackets indicate the bino- 
mial coefficient. Generally, the larger the representation 
factor is, the smaller the P-value, indicating that the 
overlap is more significant. 

Additional material 



Additional file 1: ABA-responsive genes in Arabidopsis guard cells 

This file lists the ABA-induced genes and ABA-repressed genes identified 
in Arabidopsis guard cells. 

Additional file 2: ABA-responsive genes in Arabidopsis rosette 
leaves. This file lists the ABA-induced genes and ABA-repressed genes 
identified in Arabidopsis rosette leaves. 

Additional file 3: Supplementary materials of this study. This file 
contains the supplementary materials for motif analysis and comparison 
of other ABA/hormone transcriptome studies. 

Additional file 4: New ABA-responsive genes in Arabidopsis guard 
cells. This file lists all the new ABA-responsive genes in guard cells not 
previously reported by any of the 11 previous ABA transcriptome studies 
that we evaluated. 

Additional file 5: New ABA-responsive genes in Arabidopsis rosette 
leaves. This file lists all the new ABA-responsive genes in leaves not 
previously reported by any of the 11 previous ABA transcriptome studies 
that we evaluated. 

Additional file 6: Known guard cell signaling and development 
genes from published literature. This file gives 149 guard cell signaling 
and development genes compiled from published literature which are 
classified as upregulated by ABA, down-regulated by ABA, and non- 
responsive to ABA at the transcript level. 

Additional file 7: Genes identified as ABA-regulated in guard cells 
by each of the two methods. This file lists the genes identified as ABA- 
regulated in guard cells by either the Boolean method, the linear model 
method, or both, i.e. the genes in each subsector of Figure 11 A. 

Additional file 8: Genes identified as ABA-regulated in leaves by 
each of the two methods. This file lists the genes identified as ABA- 
regulated in leaves by either the Boolean method, the linear model 
method, or both, i.e. the genes in each subsector of Figure 1 1 B. 

Additional file 9: Genes present in guard cells but not regulated by 
ABA at the transcript level. This file lists the genes that are present in 
guard cells but not regulated by ABA at the transcript level (i.e. not 
present in any of the 8 subsectors of Figure 1 1 A). 

Additional file 10: Genes present in leaves but not regulated by 
ABA at the transcript level. This file lists the genes that are present in 
leaves but not regulated by ABA at the transcript level (i.e. not present in 
any of the 8 subsectors of Figure 11B). 
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